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Method and Apparatus for Filtering and 
extending RNA alignment coverage 

5 

BACKGROUND OF THE INVENTION 
TECHNICAL FIELD 

10 The invention relates to the analysis of combination of RNA derived sequence 
fragments, hereinafter referred to as RNA. More particularly, the invention 
relates to a method and apparatus for filtering and extending RNA alignment 
coverage. 

15 DESCRIPTION OF THE PRIOR ART 

Expressed sequence tags (EST) and mRNA sequences can be aligned to 
genomic for two purposes. One purpose is to enhance gene prediction, and 
the other is to examine alternative splice forms of a gene. In the case of 
20 enhancing gene prediction, a set of non-overiapping alignments can be fed to 
a gene prediction algorithm, which uses them to give higher weight to the 
probability of there being a gene in the covered areas. With alternative gene 
splice forms, each alignment variant is fed to the gene predictor one at a time, 
which then gives greater weight to the input splice form. 

25 

These two purposes can be at odds with each other. In the case of 
enhancing gene prediction, it is desirable to extend the alignment coverage as 
much as possible to improve the gene prediction. On the other hand, each 
run of the gene prediction algorithm is CPU intensive. It is therefore 

30 preferable to run the algorithm over the same genomic region only a minimum 
number of times. In the case of alternative splice forms, running the algorithm 
repeatedly over the same genomic area is unavoidable because it can only 
take non-overiapping alignments as input. However, a problem arises when 
the different splice forms are extended and combined in such a way that a 

35 combinatorial explosion results. For example, if a gene has four areas in 
which there exist four splice variants in each area, a naive combination of 
these would result in 256 alternative splice forms. While this 
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may be conceivable, it does not seem likely. See Fig. 1 , in wliich rectangles 
represent matching blocks, i.e. exons, while lines represent gaps in the 
alignment, i.e. introns. Fig. 1 shows four aligned RNAs (A, B, C, D) from a 
gene with four areas containing four splice variations each. Permitting all 
5 possible exon combinations would yield 256 possibilities. 

These issues give rise to the need for an approach that is capable of 
extending a typically sparse RNA alignment coverage, without creating 
redundant or improbable alignments. 

10 

SUMMARY OF THE INVENTION 

The invention provides a technique that is capable of extending a typically 
sparse RNA alignment coverage, without creating redundant or improbable 

■i5 alignments. At a high level, the invention provides a two-step process. The 
first step is to combine and catenate all combinations of overlapping 
alignments that agree with each other. The second step is to extend the 

;f, boundaries of overlapping alignments that agree with their first and last exons. 

; |0 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 shows four aligned RNAs from a gene with four areas containing four 
splice variations each; 

25 Fig. 2 is a flow diagram showing a technique for filtering and extending RNA 
alignment coverage according to the invention; 

Fig. 3 shows a sample run, in which rectangles represent matching blocks, 
/.e. exons, while the lines represent gaps in the alignment, i.e. introns; 

30 

Fig. 4 shows splice variants that would be suitable for combining with other 
well-supported alternative splicings; 

Fig. 5 shows a hypothetical run of a liberal alternative splicing algorithm that 
35 creates four possible splice variants where only one is expected; 
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Fig. 6 shows a hypothetical run of a liberal alternative splicing algorithm that 
creates unsupported exon combinations; 

5 Fig, 7 is a flow diagram showing a preprocessing and filtering technique 
according to the invention; 

Fig. 8 shows a sample input; 

10 Fig. 9 is a flow diagram showing a merging technique according to the 
invention; 

Fig. 10 shows different outputs, depending on when redundant alignment B is 
;:7 discarded; 

4:;; Fig. 1 1 is a flow diagram showing an extending technique according to the 
invention; 

Fig. 12 shows a sample run of an extending step; and 

lao 

U Fig. 13 is an example showing the invention used to filter out redundant 
41 alignments. 

3 DETAILED DESCRIPTION OF THE INVENTION 

25 

Splice sites are pairs of nucleotides that indicate the location of an intron/exon 
cleavage site. The 5' site is the upstream, donor site, while the 3' site is the 
downstream, acceptor site. Together, 573' splice sites delineate the 
boundaries between exons and introns. Two alignments conflict when they 
30 disagree about the location of one or more 573' splice sites in a given 
genomic region. Two alignments agree when they overlap and there is no 
disagreement about the location of their .573' splice sites. Two alignments 
overlap if they are aligned to overlapping sections of the genomic. 
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Fig. 2 is a flow diagram showing a technique for filtering and extending RNA 
alignment coverage according to the invention. The invention herein is 
applicable to any combination of RNA derived sequence fragments, for 
example (but not for purposes of limiting the applicability of the invention to 
5 other sequences and sequence fragments) RNA derived sequence fragments 
comprising any of ESTs, partial cDNA's and full-length cDNA's. 

At a high level, the invention provides a two-step process. 

::io ■ The first step (102) is to combine and catenate all combinations of 
overlapping alignments that agree with each other. 

■ The second step (104) is to extend the boundaries of overlapping 
' alignments that agree with their first and last exons. 

! :l 

-I'r' 

Fig. 3 shows a sample run, in which rectangles represent matching blocks, 
i.e. exons, while the lines represent gaps in the alignment, i.e. introns. 

In Fig. 3, step 1 (Input), the alignment E does not merge with any other 
20 alignments. This is because its second exon does not agree with any other 
alignment. 

Step 1 serves the purpose of filtering out many redundant alignments, as well 
as combining and expanding alignments that support a self-consistent 
25 splicing. There are occasions however, when the RNA evidence for an 
alternative splicing is local to a particular part of the gene, and is not 
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able to merge with any other alignments due to its irregular set of splice sites 
(see Fig. 3, alignment E). It would be useful to have this kind of alternative 
splicing presented in the context of the entire gene, and its other splice 
variants (as in Fig. 3, alignments A-B-E and A-C-E). However, it is important 
5 to be careful not to create any exotic combinations that are not well-supported 
by the alignment evidence. 

Assessing Splice Variant Combinations 

'%0 There are many other splice variants such as shifted slice sites, and variable 
exons that it is desirable to see combined with the other variations (see 
£ Sch mucker et a/., Drosophila dscam is an axon guidance receptor exhibiting 
O extraordinary molecular diversity, Cell 1 01 , pp. 671 -684 (2000)). 

(:45 Fig. 4 shows splice variants that would be suitable for combining with other 
;2 well-supported alternative splicings in accordance with the invention disclosed 

herein. Thus, Fig. 4 shows a normal gene splicing 40, a cassette exon 41, a 

shifted splice site 42, a cassette exon with a shift 43, and variable exons 44. 

In Fig. 4, the term 'Cassette exon with shift' is explicitly listed because the 
20 inclusion of a cassette exon often causes a frame shift, leading to the 

presence of one or more nearby shifted splice sites. 

There are other splice variants however, that if combined, would lead to 
debatable alternative splicings. For example, consider an alignment that 
25 included an RNA that was only partially spliced. It would be possible to 
combine spliced RNA with it in such a way that would create many 
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undesirable splicings. 

Fig. 5 shows a hypothetical run of a liberal alternative splicing algorithm that 
creates four possible splice variants where only one is expected. In Fig. 5, 
5 there is probably only one gene splicing with all three introns. By being very 
liberal in accepting alternative splice forms, four variants can be inferred, i.e. 
the last variant (A-B-C), which is expected, plus the first three dubious 
variants (A, A-B, and A-C). 

:|0 Note that if the alignment A spanned three introns instead of two, then there 
would be seven dubious exon combinations. This situation has the possibility 
for an exponential explosion in possible splice combinations. 

With regard to the example of variable exons shown in Fig. 4, consider that 
U5 given the three alignments shown as input, it is not desirable to create 
unsupported exon combinations (see Fig. 6). As before, this situation has the 
potential for a combinational explosion in splice possibilities, most of which 
would be suspect. 

20 The solution to this problem in the herein disclosed algorithm is to preserve 
the exon configurations that are found in the input. Rather than adopt the 
notion of combining different splice variations, the idea is to extend the given 
variations. Thus, each different, shorter variation is extended in a manner 
identical to alignments A-B and A-C. This is because alignment E agrees with 

25 the others in their first and last overiapping exons. This procedure is 
discussed in greater detail below. 
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Algorithm 

Fig. 6 shows a hypothetical run of a liberal alternative splicing algorithm that 
5 creates unsupported exon combinations. The output would actually have 
seven alternative splicings, i.e. three splicings reflecting the input, and four 
fabricated splicings. 

Fig. 7 is a flow diagram showing a preprocessing and filtering technique 
^0 according to the invention. The presently preferred embodiment of the 
p invention comprises an algorithm that begins by reading a file containing 

alignments of exon annotations to genomic sequence (200) and populating an 
□ array of data structures (202), one for each alignment. Each alignment is 
Mi stored as a set of alignment blocks (204), each block representing a matching 
Ills region between the given RNA and the genomic (206). The blocks are 
Q considered exons, and the gaps between them are considered introns. For 

purposes of the discussion herein, the 573' splice sites are referred to as hard 

edges, and the two ends of the alignment are referred to as soft ends. 

20 As the alignments are read (208), some data filtering and pre-processing is 
performed. 

1 . Determine if gaps are introns or inserts (21 0). 
25 2. Eliminate single exon alignments (212). 
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3. Trim soft ends of alignments (214). 

4. Safely filter out highly similar alignments (216). 

5 First, gaps in the alignment blocks that are less than or equal to twenty 
nucleotides are considered inserts instead of introns, and the two adjacent 
blocks are subsequently combined into one. This number is a highly 
conservative, and somewhat arbitrary, educated guess at the length of the 
smallest human intron. Those skilled in the art will appreciate that other 
;|o values may be chosen in connection with practice of the invention. 

p' Second, once an entire alignment has been read and all inserts removed, the 
n alignment is discarded if it consists of only a single exon. 

115 Third, the soft ends of the alignment are trimmed by ten nucleotides. This is a 
3 heuristic for correcting a possible alignment error, in which a small number of 
terminal nucleotides are aligned to a lucky adjacent intron, instead of the 
distant correct exon. Those skilled in the art will appreciate that other values 
may be chosen in connection with practice of the invention. 

20 

Lastly, once all of the multi-exon alignments have been read and trimmed, 
there is a filtering step to cut down on running time. 

Fig. 8 shows a sample input. Multiple alignments that are very similar can 
25 lead to an exponential number of possible mergings. Filtering out similar 
alignments can be tricky. Consider Fig. 8. There are three pairs of very 
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similar alignments. In this example, alignments A and B can each merge with 
alignments C and D. The resulting alignments A-C, A-D, B-C, and B-D, can 
be merged with alignments E, F, and G to create eight combinations in a 
situation where only two combinations are necessary. A simple solution is to 
5 throw out similar alignments. Two alignments are similar if they are on the 
same strand, and have similar number of agreeing exons. This however, 
could have undesirable side-effects. 

In Fig. 8, one could eliminate shorter alignments A, C, and F because they are 
=':j[0 similar to alignments B, D, and E, respectively. This would leave a situation in 

which nothing could merge with alignment G, /.e. only A-C and B-C can align 
J; with G. On the other hand, one could eliminate alignment B, D, and E, but 
Q then nothing could merge with alignment F. 

1 115 The solution is to remove similar alignments only if their soft ends are not on 
□ opposing sides of some other alignment's hard edge. In Fig. 8, neither 
alignment C nor D would be eliminated without altering any subsequent 
merging. To keep the alignments as long as possible, whenever a similar 
alignment is eliminated, the alignment that remains is stretched to the widest 
20 possible soft end points of the two alignments. The result is that for a set of 
many similar alignments whose soft ends are not near any hard edges, only 
one alignment remains after this filtering step. Furthermore, the soft ends of 
the remaining alignment are the widest possible among all of the similar 
alignments. The last part of Fig. 8 shows the results of filtering the current 
25 example. After filtering, there are only three merging combinations (A-C-E, A- 
D-E, and A-C-G), instead of eight. 
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step 1 -Merging 

Fig. 9 is a flow diagram showing a merging technique according to the 
5 invention. The purpose of this step is to create larger alignments, while 
throwing out redundant alignments. To this end, the merge step performs a 
pairwise comparison of each overlapping RNA on the same strand (300). If 
their splice sites agree (302), then either a new, longer alignment is created 
out of the two (316; see Fig. 3 sequences A and B), or one of the alignments 

^-lO is labeled as redundant to the other (308), /.e, it is scheduled for deletion (312; 

!2 see Fig. 3, sequences C and D). 

An alignment X is redundant if it agrees with another alignment Y (304), and 
M, the soft ends of X fall completely within the soft ends of Y, inclusive (306). 
ills Redundant alignments are not deleted until after the pairwise comparisons 
3 are complete (310). 



Fig. 10 shows different outputs, depending on when redundant alignment B is 
discarded. In Fig. 10, note that a smaller alignment can be redundant with 

20 respect to another larger alignment and still be able to merge with other 
alignments that the larger alignment cannot. Redundant pieces are still 
available for the pairwise comparisons, even though they are redundant with 
other alignments. This is why performing all pairwise comparisons is 
important. The significance of noting redundant alignments is that they should 

25 eventually be discarded because they are redundant. 
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Notice in Fig. 10 that if alignment B were completely ignored because it is 
redundant with alignment A, then it would not be available to merge with C. 

Each newly created, merged alignment is also checked against the other, 
5 unmerged alignments in the same way (318): either merging once again, 
becoming labeled as redundant, or causing other alignments to be labeled 
redundant. Whenever an alignment, merged or not, has been compared to all 
other alignments (320), and is neither redundant nor merged it is placed on a 
list of Done alignments (314). Once all comparisons have been finished 

10 (322), the Done alignments are again compared pairwise to check for 

';f; redundancies (324, 326, 328). 

.C' The possibility of redundant alignments within the Done list is seen in Fig. 4. 
. There, alignment A merges with C, and the combination can not merge any 

further, so it is done. Alignment A also merges with D, and that combination 
;5 is also done. Now both the A-C and A-D combination are on the Done list, yet 

the A-D combination is redundant. 

In Table 1 below, the algorithm is guaranteed to terminate because the items 
20 on the todo list always begin their comparisons in the cdnaArray after the 
index where the merging took place. Therefore, attempting to compare the 
same merged alignment ad infinitum is impossible. 

The merge portion is complete when the Done list has been cleared of 
25 redundant, merged alignments. 
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Table 1. Pseudo-code for merging algorithm 



/ * cdnaAray = array of cDNA alignments stored by increasing 

* start points 

5 * todo list = a list of merged alignments that are 

scheduled to be compared to the other * 
alignments in cdnaArray. Associated with each 

* merged alignment is an integer 

* array index indicating where in cdnaArray the 
10 * comparisons should begin. 

V 

bool merged; 
,,,,, for I = 0 to sizeOf (cdnaArray) { 
!|| for j = (I + 1 ) to sizeOf (cdnaArray) { 

lis merged = 0; 

%j if ( overlap (cdnaArray[l], cdnaArrayU])) { 

F merged I = merge (cdnaArray[l], cdnaArray[j], todojist); 

a } 

:Jo if (Imerged) { 

m push(DoneJist, cdnaArray[l]); 

'Pi else {/* there must be a mered alignment on the todojist V 

while (todo_list !=NULL) { 
25 todo = pop(todo_list); 

merged = 0; 

for k = todo,dtart to size Of (cdnaArray) { 
if (overlap (todo, cdnaArray [k])) { 

merged l=merge (todo, cdnaArray [k], todojist); 

30 } 
} 

if (Imerged) { 
push(Done-list, todo); 

} 

35 } 
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step 2-Extending 

5 Fig. 1 1 is a flow diagram showing an extending technique according to the 
invention. Fig. 12 is a sample run of Step 2, extending. Combination (A-D) is 
possible only with the presence of alignment (D). It cannot be inferred from 
the other alignments. 

^'|o This step also involves a pairwise comparison of all the remaining alignments 
p (400). For each pair, the left-most and right-most overlapping exons are 
g considered (402). If there are no conflicts on any of their hard edges (406), 

r% then the short ended alignments are extended in such a way that they match 

u the longer alignments (41 2). 

m 

C3 Notice in Fig. 12, after step 2, there would be no alignment that represents the 
exon combination of alignments B and C. This is because alignment C's final 
exon conflicts with B, and therefore the algorithm does not use this 
combination. If there were also an alignment D (shown in parentheses in Fig. 
20 12) that was similar to C, only its final exon agreed with alignment B, then all 
four exon combinations would be considered. 

The motivation for considering all pairs of alignments is that once an 
alignment has been extended, it has new opportunities to extend, and be 
25 extended by, other alignments. If an alignment can be extended, then a new 
alignment is created (414), and it is further compared against all 
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overlapping alignments. If an alignment cannot be extended, it gets placed on 
a Done list (408), as before. After all comparisons are complete, the Done list 
is again purged of redundant alignments (410). The pseudo-code for this step 
is almost identical to the pseudo-code for the merging step, and is therefore 
5 omitted. 

What remains after Step 2 is a minimal set of splice variant alignments, 
extended as far as they can given the alignment data. 

The presently preferred embodiment of the invention thus comprises a 
J2 conservative, yet effective technique that addresses the issue of combining 
5; and expanding aligned RNA's to create a concise set of inputs for a gene 
Q prediction algorithm. The preferred embodiment of the invention is capable of 
u extending small RNA alignments that represent an alternative splicing into 
1115 much larger alignments, consistent with the other well-supported splicings. 
Q This is done in a manner that does not create new exon combinations that are 

not supported by the given data. The preferred embodiment of the invention 

also outputs a set of alignments, each of which represents a different set of 

splice site combinations. 

20 

Although the invention is described herein with reference to the preferred 
embodiment, one skilled in the art will readily appreciate that other 
applications may be substituted for those set forth herein without departing 
from the spirit and scope of the present invention. In addition to the 
25 algorithm's ability to derive alternate splice form combinations, it has the 
capacity to filter out redundant alignments. As more and more RNA fragments 
are sequenced and recorded in genomic databases, the need for automated 
filtering of this data becomes greater and greater. Thus, in a situation where 
hundreds of RNA alignments overlap within the same genomic 
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region, the algorithm could be used to cull out the possibly few useful splice 
forms, while purging the redundant alignment data (see Fig. 13). Accordingly, 
the invention should only be limited by the Claims included below. 
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